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ABSTRACT 

Lensing measurements of the shapes of dark matter halos can provide tests of gravity 
theories and possible dark matter interactions. We measure the quadrupole weak lens¬ 
ing signal from the elliptical halos of 70,000 SDSS Luminous Red Galaxies. We use a 
new estimator that nulls the spherical halo lensing signal, isolating the shear due to 
anisotropy in the dark matter distribution. One of the two Cartesian components of 
our estimator is insensitive to the primary systematic, a spurious alignment of lens 
and source ellipticities, allowing us to make robust measurements of halo ellipticity. 
Our best-fit value for the ellipticity of the surface mass density is 0.24 db 0.06, which 
translates to an axis ratio of 0.78. We rule out the hypothesis of no ellipticity at the 4cr 
confidence level, and ellipticity < 0.12 (axis ratio > 0.89) at the 2cr level. We discuss 
how our measurements of halo ellipticity are revised to higher values using estimates 
of the misalignment of mass and light from simulations. Finally, we apply the same 
techniques to a smaller sample of redMaPPer galaxy clusters and obtain a 3cr mea¬ 
surement of cluster ellipticity. We discuss how the improved signal to noise properties 
of our estimator can enable studies of halo shapes for different galaxy populations 
with upcoming surveys. 

Key words: cosmology: observations - dark matter; gravitational lensing: weak; 
galaxies: haloes 


1 INTRODUCTION 

While isotropic dark matter distributions can accurately de¬ 
scribe galaxy-galaxy lensing statistics of halos (Mandelbaum 


let al. 2005 Johnston et al. 2007 Sheldon et al. 2009b), 


there is a theoretical expectation for individual halos to have 
significant ellipticity (Jing & S uto| |2002). A number of au¬ 


thors have attempted to measure this higher-order signal on 


galaxy-scale i 

jHoekstra et al. 2004 Mandelbaum et al.|2006 

Parker et al. 

2007 van Uitert et al. 2012) and cluster-scale 

halos (Evans & Bridle 2009 Oguri et al. 2010). The cluster 


measurements, where the intrinsic signal is expected to be 
much stronger, have been more successful. Using 4,300 clus¬ 


ter lenses Evans & Bridle (2009) obtain a 2.7cr detection of 


signal. Furthermore, using a subset of 18 clusters that are 
well fit by an elliptical NFW model, Oguri et al. (2010) esti¬ 


mate a 7cr detection of ellipticity for those specific clusters. 

Missing from the literature until now is an attempt to 
detect and characterize the ellipticity of group-size halos, 
such as those in which Luminous Red Galaxies (LRGs) re¬ 
side. With intermediate size between galaxies and clusters, 
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these LRG halos have larger intrinsic lensing signal than 
the galaxy samples and are larger in number than the clus¬ 
ter samples. However, as with the galaxy halos, these have 
neither enough member galaxies nor enough unstacked lens¬ 
ing signal to make an initial guess for the halo major axis. 
We must rely on the ellipticity of the light distribution in 
order to perform an aligned, stacked lensing measurement. 
This approach has both advantages and disadvantages: on 
the positive side our method has promise in constraining 
Modified Newtonian Dynamics ( Milgrom|1983 ) and related 
modified gravity theories ( Khoury[2 015), which rely on the 
distribution of baryonic matter alone to produce every ob¬ 
servable. Such theories will have difficulty producing suffi¬ 
cient anisotropic lensing signal at scales much greater than 
the halo’s member galaxies. On the other hand, this method 
is susceptible to a possible systematic alignment of lens and 
source ellipticities (Mandel baum et al.|2006 ) and constraints 
must take account of possible misalignment of the dark mat¬ 
ter and light distributions ( Adhikari et al.||2015 ). Nonethe¬ 
less, we show that our method and resulting lower bounds 
on halo ellipticity are robust even in the presence of these 
systematics. 

The measurement in this paper is along the lines of our 
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2 Clampitt and Jain 


previous work producing lensing detections that go beyond 
the description of isotropic dark matter halos to measure 
higher order structures of the cosmic mass distribution. In 


Clampitt et al. (2016) we detected the weak lensing signal of 


filaments, and in Clampitt & Jain (2015) we constrained the 
density profile of cosmic voids. For all these measurements, 
filament lensing, void lensing, and now elliptical halo lens¬ 
ing, we have relied heavily on the Sloan Digitial Sky Surve}Q 
(SDSS) shear catalogs of Sheldon et al. (2009a). While new 
lensing measurements from the Dark Energy Surve}p](DES) 
are starting to be released (Melchior et al. 2015; Cha ng et| 


2015; Vi kram et al.]|2 
rveyj and the Subaruj 


al. 2015) Vikram et al.])2015) (as well as the Kilo Degree 


Survey and the SubarijjHSC survey), there is still more to 
be gleaned from the previous generation data of SDSS. Fur¬ 
thermore, the measurement techniques developed in these 
works will be useful to apply to DES and other cutting edge 
datasets in the near future. 

In § 2 we describe our model for the shear from ellip¬ 
tical halos as well as potential systematic errors. In § 3 we 
describe our data for lens and source galaxies and also spec¬ 
ify our observables. Then in § 4 we show the results for the 
measurement and constraints on halo ellipticity, compare to 
previous estimators used in the literature, and repeat the 
measurement with redMaPPer clusters. In § 5 we discuss 
systematic shear tests and the impact of galaxy-halo mis¬ 
alignment. Finally, in § 6 we summarize our main conclu¬ 
sions, compare to previous work, and discuss implications 
for theories of gravity and dark matter. 


2 MODEL 


2.1 Elliptical halo shear signals 


We use the results of Adhikari et al. (2015), whose model 
for the effect of halo ellipticity on weak lensing observables 
will be useful in interpreting our measurements. This model 
for the surface density of elliptical halos uses a multipole 
expansion: 

S(i?, 9) oc R vo [1 + e(— ?7o/2) cos 29 + 0(e 2 )] (1) 

= £ 0 (-R) + E 2 (-R) cos 2 # + ... ( 2 ) 


where 770 (R) = d log Eo /d log R < 0 and we assume the coef- 
fecient of the quadrupole —erio/2 <C 1 , justifying the neglect 
of higher orders in the expansion. Here R is the projected 
distance from the center of the halo, and 0 is the angle rel¬ 
ative to the halo’s major axis. We set 


—2 E 2 (R) 

rio(R)Z 0 (R) ’ 


(3) 


thus allowing the quadrupole E 2 to be completely deter¬ 
mined by the monopole Eo, up to a proportionality factor 
e, the magnitude of the ellipticity. This ellipticity is related 
to the axis ratio of the projected mass distribution, g, by 
e = (1 — q 2 )/(l + q 2 ). Note that our ellipticity definition is 
bounded on the interval [ 0 , 1 ] and thus differs by a factor of 
2 from Adhikari et al.| (2015). 


1 http://www.sdss.org 

2 http://www.darkenergysurvey.org 

3 http://kids.strw.leidenuniv.nl 

4 http://www.subarutelescope.org/index.html 
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Figure 1. The quadrupole shear pattern produced by an elliptical 
dark matter halo. The x-axis of the Cartesian coordinate system 
is aligned with the major axis of the galaxy light distribution, 
assumed to be aligned with the major axis of the dark matter 
halo. Our sign convention for the Cartesian shear components is 
shown at the right. Note that the monopole shear (which is purely 
tangential) is not contributing to the pictured shear pattern. 


Adhikari et al. (2015) shows that the tangential and 


cross-components of shear from the quadrupole are given by 


E C rit 7+ 
E C rit 7x 
where 


= (e/ 2 ) [Eo (R)rio(R) - Ii(R) ~ h(R)} cos 2 0 (4) 

= (e/ 2 ) [-h (R) + h (R)] sin 20 , (5) 


h(R) - w 


( 6 ) 

(7) 


J -R ,3 £o(-R , )?7o(-R , )di^ , , 

hr) = j f* y) drf. 

(The lensing weight function E cr i t is defined later in Eq.[2l)) 
Note that throughout this paper we are concerned only with 
shear from the quadrupole: all shear variables 7 +, 7 x, etc., 
refer to shear from the quadrupole and have no contribution 
from the monopole. 

In § |4.3| we take a careful look at weighted estimators of 
tangential and cross shear that have been studied elsewhere 
in the literature (Natarajan V Refregier 2000; Mandelbaum 


et al.|2006 van Uitert et al.|2012 ). However, for most of this 

work we use Cartesian estimators, measured with respect to 
a coordinate system with the x-axis aligned with the lens 
LRG’s major axis. Our sign convention for these Cartesian 
71 and 72 components is shown in Fig. [l] In this coordi¬ 
nate system, the two shear components are related to the 
tangential and cross-shear by 


71 CM) 

Ecrit ^71 (R,0) 


and 


72 (R,0) 

Ecrit 72 (R,0) 


— 7 + (R, 0) cos 26 + 7 x (R, 0) sin 20 ( 8 ) 

(e/4) [(2/i (R) - Eo(RMR))cos4<9 + 
2/ 2 (R) — E 0 (R)?7o(R)] , (9) 

— 7 +(R, 0) sin 20 — 7 X (R, 0) cos 20 ( 10 ) 
(e/4) [2/i (R) - E 0 (R)? 7 o(R)] sin 40 . (11) 


Studying these equations, we see that the angular depen¬ 
dence of the shear goes as cos 40 or sin 40. Thus there is 
a sign change in both components after every angle 7 r/ 4 , 
and moving around the circle the shear signal from elliptical 
halos transitions between regions where 71 then 72 alter¬ 
nately dominate. The resulting shear pattern of Eqs. (|9 11) 
is sketched in Fig. [T| 
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Halo Ellipticity 3 


There are two minor exceptions to the shear direction 
of Fig. [l] both are a result of the two terms without 0 de¬ 
pendence in Eq. for 71 . First, for a small part of our 
measurement range (below R ~ 0.14 Mpc /h, see Fig. [ 2 ]) the 
sign of 71 is negative for all 6. Second, the zero-crossing of 
71 for different 0 depends somewhat on scale R. This leads 
to a small mixing of positive and negative shears in the sep¬ 
arate ranges of 6 we select below. However, we have checked 
that this only leads to an appreciable loss of S/N in a small 
range around R ~ 0.14 Mpc /h: at scales R < 0.1 Mpc /h and 
R > 0.3 Mpc /h the loss of S/N due to mixing positive and 
negative shears in one bin is less than 10%. Both of these 
effects are taken into account in our modeling so that our 
constraints on halo ellipticity are not biased. Note that while 
we are only concerned with the total signal from the halo 
quadrupole, the origin of the constant shear term has been 
studied by Bernstein & Nakajima (2009) who determined 
that it comes from quadrupoles external to the source radius 
R. (Note that at very small radii R , most of our elliptical 
halo is an “external quadrupole,” since this mass is outside 
the radius of interest.) This distinction between terms com¬ 
ing from quadrupoles internal and external to the source 
radius may be useful in future higher signal-to-noise mea¬ 
surements, e.g., in discriminating between the quadrupole 
of the dark matter halo and that of correlated large-scale- 
structure. Since our measurements focus on small scales, in 
this work we do not include a quadrupole from correlated 
large-scale-structure in the model. 

These sign changes suggest grouping together the 71 
component of background sources in the regions where 
cos(40) has the same sign, i.e., 

rn/8 

/ E cr it 7 i(R, 6) dO and 

J - ?r/8 

/»37 t /8 

/ E crit 71 (.R, 0) d6 , 

n J n/8 

where we have introduced the notation of AE in analogy to 
tangential shear, for which AE = E cr it 7 +- We do the same 
for the Cartesian cross-component, 72, except here we group 
together same-sign regions of sin(4$): 

- f ' S cri t72(#,<?)d0, (14) 

k Jo 

- I"' £ crit 72(#,0)d6>. (15) 

n Jir/4 

In order to capture the information on halo ellipticity in the 
shear field, we measure these four functions from the data. 
Note that the above integrals only cover 1/4 of the circle 
for each shear component, but including the other regions 
offset by an angle 7r/2 gives equivalent results, as can be 
seen in Fig. [T] When measuring these functions, in order to 
increase statistics we exploit these symmetries and average 
together all regions with the same sign of the shear (see the 
estimators in § [ 3 ]) . 

These estimators also have the useful property that 
purely tangential shears (the lensing monopole) are nulled 
automatically. This can be seen from Fig. [l] consider two 
source galaxies located at (x,y) positions (R, 0 ) and ( 0 , R ) , 
i.e., in the regions where we measure AE^ as in Eq. (12). 
Both have monopole contributions of 7 * > 0 which are 
purely tangential and equal in magnitude (since they are lo- 


AS [~\R) = - 


AE [ +) (R) = 


( 12 ) 


(13) 


AE ( ~\R) = 


AE< +) (fl) = 


cated at the same distance R from the LRG center). Trans¬ 
forming to Cartesian shears ( 71 , 72 ), the two sources have 
(— 71 , 0) and ( 71 , 0), respectively. When averaged in the same 
bin, the net monopole contribution from these two sources is 
zero. More generally, whenever two sources are offset by 90°, 
relative to the LRG center, the Cartesian components of the 
monopole shear cancel. Assuming background sources are 
distributed isotropically around the lens, each source used 
in any of the four estimators in Eqs. HU will always have 
a paired source offset by 90°. (Note that we have checked 
that our sources are indeed distributed isotropically.) For 
more detailed discussion of this nulling of monopole shear, 
and its usefulness in constructing an estimator of filament 


lensing, see Clampitt et al. (2016). 


Note that in the absence of systematics, the form of 
Eqs. ([9 


11 ) suggests estimators weighted by sin 4 0 and 


cos A0. The weighting is not optimal since it actually loses 
signal in 71 by averaging out the terms in Eq. that have 
no angular dependence. In addition, the presence of signifi¬ 
cant galaxy-halo misalignment (§ 5.2) introduces an error in 


the choice of the angle 9 for each lens-source pair and dilutes 
the stacked signal. We find that for our sample the simple 
sin 4 6 and cos 4 0 weighting does not lead to an improvement 


in the signal-to-noise (see § 5.3). We leave for future work 
the exploration of estimators that may improve on the ones 
presented here. 

Our model predicts the shear quadrupole from the 
monopole, up to an overall factor of the ellipticity. Thus 
we use the parameters of these LRG halos which have been 


determined by Mandelbaum et al. (2008a) using isotropic 


galaxy-galaxy lensing: the average mass and concentration 
are M — 4 X 10 13 M®//i and c = 5. These parameters, along 
with our assumed NFW halo profile ( Navarro et al. |[l997 ) 
then fully determine the signals of monopole E o(R) and 
quadrupole, Eqs. (12)-(15). We discuss below the impact 


of varying M and c on our ellipticity estimate. 


2.2 Systematic shear 


For the SDSS survey, Mandelbaum et al. (2006) found that 


nearby galaxy shapes tend to be slightly aligned with each 
other. There are several factors contributing to this align¬ 
ment. First, the PSF has a net large-scale alignment with 
the survey scan direction, and any incomplete correction of 
this PSF will tend to affect nearby galaxies in the same way. 


In addition, Huff et al. (2014) found that specific declina¬ 


tion ranges have systematic shear due to an incorrect PSF 
model in the r-band camcol 2 CCD; see especially the dis¬ 


cussion around their Fig. 8 . (Note also Sec. 4.4.3 of Reyes et 


al. ( 2012 ), which studied the effect of this systematic on the 


shear catalog of Mandelbaum et al.|(2006) in more detail. 


Both effects imprint a slight alignment between our lens and 
source ellipticities. For isotropic galaxy-galaxy lensing, these 
systematics are not a worry because they are largely can¬ 
celled out along with other additive errors. 

One significant advantage of using the Cartesian shear 
components described above to measure halo ellipticity is 
that this major systematic error only affects the AE^ +//_ ^ 
functions, leaving AE^ ^ untouched. Once one averages 
over many lens-source pairs, the only preferred orientation 
is the major/minor axis of the lens which are aligned with 
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4 Clampitt and Jain 


the x- and y-axes. A systematic error can contribute to 
the estimated shear along these axes, and thus contami¬ 
nate the 1-component, as discussed below for SDSS. How¬ 
ever by symmetry, the net contribution of a systematic to 
the 2-component of the shear must vanish since both signs 
are equally likely, unlike the quadrupole signal. Previous at¬ 
tempts at measuring halo ellipticity used the tangential and 
cross-components of the shear, which can both be contam- 
inted by systematic errors. One way to see this is the 1- 
component has non-zero projections along both the tangen¬ 
tial and cross-components (except at special points such as 
0 = 0). Note that in addition to these arguments from sym¬ 
metry, in Fig.|7|we have tested independently that ^ 

has no contribution from this systematic. 

The elliptical halo shear measurement of Mandelbaum 
et al. ( 2006| ) is dominated by this systematic for at least one 
lens sample (the red “L4” lenses in Fig. 7 of that work). For 
this sample, the systematic has a roughly 1/R power law 
dependence. Thus we parameterize the systematic shear as 


et al. (2010) for more details on the LRG catalog, Sheldon et 


A£ syst oc ^ 

R 


(16) 


and marginalize over the amplitude, which is required to be 
positive according to our shear convention. This can be seen 
from Fig. [I] where an alignment of the lens galaxy and source 
whiskers will contribute positively to 71. We have also tried 
other power law models, such as 1 /R 2 , for the systematic 
shear and find our results are not sensitive to the power 
chosen. Our full model for the A£i components is therefore 
A£^ +//_ ) + A£ syst . Note that while our main results involve 
modeling and marginalizing over the systematic, in § |5.1| 
and Fig. [7] we compare with model-independent tests of the 
systematic shear and find similar results. 

Finally, we comment briefly on the possibility of a sys¬ 
tematic anti-alignment of lenses and sources (negative am¬ 
plitude for Eq. (16)). The incorrect PSF model in the r- 
band camcol 2 CCD is a large effect that always acts to 
align nearby galaxies. For the other camcols, since lenses 
and sources have different magnitudes and sizes, one could 
imagine the possibility that either the lens or source would 
be overcorrected and the other undercorrected. Thus, when 
performing the fits we will also test the results when drop¬ 
ping the requirement that the systematic amplitude be pos¬ 
itive. 


3 DATA AND MEASUREMENT METHOD 

We use shear catalogs and photometric redshift probability 
distributions that are almost identical to those in ISheldonl 


et al. (2009a) and Sheldon et al. (2012), respectively. The 


catalog contains 34.5 million sources, with a redshift distri¬ 
bution peaking at z ~ 0.35 and a long tail towards higher 
redshifts. In order to find the major axes of the LRG light 
distributions, we match the shear catalog of Sheld on" et al. | 
(2009a) to the DR7-Full LRG catalog of Kazin et al.| ( |2010| k 
This is necessary to align our Cartesian coordinate systems 
for each lens prior to performing a stacked lensing measure¬ 
ment. This results in 69,286 matching galaxies with redshifts 
0.16 < z < 0.47 in our lens sample. Essentially all of these 
matched LRGs have quality shear measurements, based on 
the small shape measurement error of each galaxy. See [Kazin| 


al. (2009a) for more details on the shear catalog, and Shel¬ 


don et al. (2004) for more detailed descriptions of the shear 


measurement method. Note that although the discussion of 
systematic shear in Sec. |2.2| focuses on the shear catalog of 
Reyes et al.| ( 2012| , the same PSF models were used for our 
shear catalog (Erin Sheldon and Rachel Mandelbaum, pri¬ 
vate communication). Since the source of systematics was 


primarily the PSF models (Huff et al. 2014 Reyes et al. 
2012), the results are also relevant for our work. 


To perform the measurement, we begin by rotating the 
RA and DEC coordinates of all nearby background galaxies 
into a Cartesian coordinate system with x-axis along the 
LRG major axis (as pictured in Fig. [I]). The source catalog 
has shear components (71,72) which are also defined with 
respect to the RA, DEC coordinate system. We simultaneously 
rotate these into (71,72) components defined with respect to 
the ( x , y ) system, keeping in mind that shears rotate twice 
as fast as galaxy positions (7 = yT 2 * 0 for a rotation by O'). 
Our lensing observable is then given by 


y. 

__ r 8 \ ^3 

A E* (R) = - 


(<e c 7>.(zl)) \ ar ) 


E, 


(17) 


where the sum goes over all sources j that fall in the specified 
region, different for each shear component k and sign s : 

k = 1, s = T* 0 < (6j + 7r/8) < 7r/4 
k = 1, s = — 1 7 t/4 ^ (0j + 7r/8) < 7r/2 

k — 2, s — +; 0 ^ 0j < 7r/4 

k = 2, s = tt/4 ^ 0j < tt/2 (18) 

and where each estimator also includes the sources in sym¬ 
metrical regions shifted by 7r/2, 7r, and 37r/2, as shown in 
Fig. [l] Note that the Cartesian shear of the j’th source 
galaxy, jk,j, should be distinguished from the predicted 
shears 71 and 72 referred to in our model Eqs. ( 4|15 ). The 
summation y. runs over all the background galaxies in the 
radial bin R, around all the LRG positions. The weight for 
the j -th galaxy is given by 


[< 


(2l) 


Tshape ^meas ,j 


(19) 


The method for these weights follows that used for tangen¬ 
tial shear measurements ( Mandelbaum et al.|2 008b; Sheldon 


et al.||2012[ Mandelbaum et al.||2013 ). We use cr s h ap e = 0.32 

for the rms intrinsic ellipticity and cr mea s,j denotes measure¬ 
ment noise on each background galaxy ellipticity. (E~[ t )^. is 
the lensing critical density for the j- th source galaxy, com¬ 
puted by taking into account the photometric redshift un¬ 
certainty: 


( S crit),- (*l) = 


f 


d*Ej t (*L, 


s)Pj( z s 


( 20 ) 


where zl is the redshift of the LRG and Pj(z a ) is the prob¬ 
ability distribution of the photometric redshift for the j-th 
galaxy. Note that £~J t (^L,^ s ) is computed as a function of 
lens and source redshifts for the assumed cosmology as 


£crit (^L 5 ^s) — 


Da(z s 


47 iG Da(^l)Da(^l, £ s 


( 21 ) 
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Figure 2. (Left panel): The measured values of AS, } (black points) and best-fit elliptical halo model (solid line), along with the same 
for ASj f (green triangles and dashed line). (Right panel): Same as the left panel, but showing the data and predictions for AE, 2 • For 
all cases, the model is able to describe the data very well. As predicted from Eqs. § and GD> the constant (no angular variation) shear 
term causes AE^ ^ to have the largest signal. The best fit value of ellipticity e is 0.24, while the systematic shear in our model has a 
best-fit value which is close to zero for all scales R on this plot. 


and we set = 0 for z s < in the computation. 

The angular diameter distances from the observer to source, 
observer to lens, and lens to source are given by Da(z b ), 
Da(zl), and Da(zl,z s ), respectively. 

The error bars and covariance are obtained with the 
jackknife method of |Norberg et al. (2009), using a set of 134 
jackknife patches. Note that given the size of the error bars 
in the following section, a few percent uncertainty in the 
multiplicative bias is not significant. In addition, the effects 
of additive bias will be treated in detail in § |5.1| Thus we 
directly use the jackknife covariances in all plots and fits. 


models will be needed with higher signal-to-noise measure¬ 
ments from upcoming survey data. 

The other components are predicted to have lesser over¬ 
all signal-to-noise (S/N) than AE^, and this is consistent 
with what we see in the data. Since the noise in all compo¬ 
nents are expected to be dominated by shape noise at these 
scales, and each has the same shape noise, the components 
with lesser predicted signal should also have lesser S/N. The 
predicted AE^and AE^ have identical amplitude and 
lower S/N, and AE^ has the least. 


4 RESULTS 

4.1 Measurement with Cartesian components 

The results of the measurement and best-fit model predic¬ 
tion are shown in Fig. [5] At R ~ 0.3 Mpc /h the magnitude 
of the measured signal is ~ 5M©/i/pc * 2 for AE^ - \ Compar¬ 
ing to measurements of the monopole signal for this LRG 
sample (Mandelbaum et al. 2008a, 2013), the magnitude of 

each component of the quadrupole is ~ 10 times smaller. By 

2 Mpc /h the measured signal is consistent with zero. Since 
this is the transition s cale between the 1- and 2-halo regimes 
(Velander et al. 2014) for halos of mass ~ 4x 10 13 M©//i, our 


4.2 Constraints on halo ellipticity 

In Fig. [3] we show our constraints on the halo ellipticity and 
the systematic alignment amplitude. The 2D constraints in 
the top panel show the systematic amplitude (see Eq. [Tg] ) is 
very small for our LRG lens sample, with a best-fit ampli¬ 
tude less than ~ 0.1 x 10 12 M©/Mpc. Thus, at R = 1 Mpc /h, 
the systematic has magnitude less than 0.1M Q /i/pc 2 in 
units of AE, and at all scales it is much smaller than our 
modeled elliptical halo signal. Marginalizing over the am¬ 
plitude of the systematic, we rule out the hypothesis of no 
ellipticity e at 4a, and e < 0.12 at 2a. These constraints can 
be read off from the lower panel, which shows the cumulative 
probability in the low e tail of the posterior. Possible mis¬ 
alignment between the galaxy and dark matter distributions 


simple 1-halo modeling of Eo(-R) is sufficient to describe the 
signal for the purposes of this study. However, more detailed 


only dilutes our measurement (see the discussion in § 5.2), so 
these lower bounds on halo ellipticity still hold in the pres- 
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Halo ellipticity, e 



Halo ellipticity, e 


Figure 3. (Top panel): Combined constraints from all four func¬ 
tions A '5^12 ^ on the halo ellipticity e and systematic shear am¬ 

plitude. For allowed values within the 2a contours, the magnitude 
of the systematic is at most 1% of the signal over all measured 
scales. (Middle panel): Probability of the halo ellipticity given the 
data, marginalized over the systematic amplitude. Constraints are 
pictured for our entire data vector (solid line) and only the compo¬ 
nents that are not dependent on the systematic model, AE^^ ^ 
(dashed line). (Lower panel): Cumulative probability distribution 
of the middle panel, corresponding to the p-value of a given ellip¬ 
ticity e. The l-4cr confidence levels (horizontal dotted lines) show 
that we rule out the hypothesis of no halo ellipticity e at 4a, and 
e < 0.12 at 2cr. These lower bounds hold even in the presence of 
misalignment between the light and dark matter profiles. 



Figure 4. The tangential- (blue points) and cross-shear (red 
triangles) estimators used in [Mandelba um et al.| ( f2006| >, as well 
as the corresonding best-fit model prediction to each (solid and 
dashed lines, respectively). These also give a significant signal 
with our lens and source sample, but subtraction of the two com¬ 
ponents to deal with systematic shear will also remove real signal 
from the elliptical halo. 


ence of misalignment. The best-fit model has e = 0.24±0.06. 
(Note that 0.24 refers to the mean of the posterior and 0.06 is 
the lcr confidence interval.) With a reduced % 2 of 25.1/22, 
the e = 0.24 model is a good fit to the data. In contrast, 
taking the null hypothesis model of no halo ellipticity the 
reduced x 2 is 41.7/24, a very poor fit. It should be noted 
that this best-fit value may be an underestimate due to mis¬ 
alignment of the light and dark matter position angles. We 
discuss this possibility in detail in § |5.2| Note that the covari¬ 
ance between the 24 data points in Fig. [2] is small but non- 
negligible: thus the full covariance has been used in these 
constraints. 

Also shown in Fig.jsjis the constraint from the AE^ +// ^ 
observable alone, i.e., not including the higher signal-to- 
noise AE^ + ^ - ^ constraints. Even if our model for the sys¬ 
tematic is somehow flawed, this observable is unaffected by 
the systematic and it gives a 2cr constraint by itself. In 
§ |2.2| we discussed the motivation for not allowing a sys¬ 
tematic anti-alignment (which will be somewhat degenerate 
with the constant-with-angle terms of Eq. @) in our fits. 
Nonetheless, if we redo the fits allowing a negative ampli¬ 
tude of Eq. (16), the mean of the systematic amplitude is 
then ~ (—0.35 ± 0.4) x 10 12 M Q /Mpc, still consistent with 
zero. Furthermore, given our tests in § |5.1| with a foreground 
“source” sample that has no lensing contribution, we think it 
unlikely that the systematic produces a net anti-alignment. 
In the following section, we repeat the measurement using 
estimators previously used in the literature, and make com¬ 
parison to previous results of Mandelbaum et al. (2006). 
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Figure 5. Same as Fig. [2] but for redMaPPer clusters. The measurement is still signal dominated, but for these clusters the systematic 
(dotted line) is not negligible with an amplitude ~ 20 — 40% of the halo ellipticity signal (dash-dot line). The total signal (solid and 
dashed lines) is 2-3 times larger than for the LRG sample. 


4.3 Measurement with tangential components 

The following estimators were developed by 


Natarajan &; 


Refregier (2000) and applied to survey data by 


et al. 


Mandelbaum 


(2006): 


/AE(fl) = 


£^(<2 crit } j (zl)) 7 +,j cos 

2 V. Wj cos 2 20j 


20n 


( 22 ) 


for the tangential component of shear and 


/ 45 AE(Ji) = 


£j w o (( s crit)j ( z l)) 7X,i sin 29j 

2 Y j w j sin 2 2 0j 


(23) 


for the cross-component. In this section we repeat the mea¬ 
surements on our lens and source samples using these es¬ 
timators, and compare the results with Mandelba um et al.| 
(2006). Given our signal in Fig. [2] we would expect both 
estimators to give positive signals when repeated using our 
lens and source samples. Indeed, this is what we see in Fig. [4] 
for /AE, while the results for AE are more noisy. The 
positive sign is partly a matter of sign conventions: ours are 
given by Fig.[l]and Eqs. ([8][ToJ - Note that our Eqs. (22) and 
(23) include an extra factor of 2 in the denominator relative 
to Eq. Pi in order to match the convention of Mandel¬ 


baum et al. (2006). In Fig. [4] we plot 2/AE and 2/45 AE so 


that the signal amplitude can be compared fairly with our 
fiducial measurements in Fig. [2] Mandelbaum et al. (2006) 
also included a factor of the lens ellipticity in the estimators 


Eqs. (22 23); we have found that including this factor has 


a negligible effect on the measurement. We have chosen to 


leave it out for simplicity and to better match the form of 
our own estimator equations. 

The most comparable lens sample of Mandelbaum et al. 

(2006) is the red “L6” lenses, for which they find an ~ 2 — 30- 
hint of alignment between the elliptical halo and lens galaxy 
light distribution. Our stronger detection is likely due to a 
few effects. While their sample has ~ 50% more galaxies 
than our LRG sample (~ 107,000 compared to ~ 70,000), 
the lack of spectroscopic redshifts likely dilutes their signal 
somewhat. Perhaps more importantly, our estimator does 
not sacrifice signal in order to remove systematics. |Mandel-1 


baum et al. (2006) argue that both estimators have an equal 


contribution from the systematic shear, and so they use the 
difference /AE — /4sAE to produce a halo ellipticity mea¬ 
surement that is free of this systematic. They also note that 
this subtraction is not ideal, however, as there is also signal 
in /45AE (see Fig. 2 of that work). When subtracting the 
two estimators, the systematic will be cancelled, but some 
real signal will also be removed at the same time. This can 
be seen in Fig. [4] where above R ~ 0.1 Mpc /h both estima¬ 
tors have positive signal from the elliptical halo. At scales 
~ 0.5 Mpc /h and larger, the two have similar magnitudes as 
well, so that subtraction will remove most of the real signal 
as well as the systematic shear. Thus we continue to focus 
on the Cartesian shear estimators of halo ellipticity in this 
paper. 


4.4 Measurement with redMaPPer Clusters 

As a check of our measurement method and interpretation, 
we attempt the same measurement with redMaPPer clusters 
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8 Clampitt and Jain 


(Rykoff et al. 2014[ Rozo et al. 2015). We consider only 


those clusters where the most-probable central galaxy has a 
high probability of being the true center (P cen > 90%). In 
addition, we use only clusters with redshift 0.1 < 2 : < 0.3, 
removing high-redshift clusters with noisier detections. We 
then match these galaxies to our shear catalog, yielding the 
direction of the galaxy major axis. Of the original ~ 26,000 
redMaPPer clusters, only ~ 2,700 remain after these three 
cuts. 

As for the LRG sample, we assume the cluster den¬ 
sity follows an NFW profile and the quadrupole of the sur¬ 
face density is completely determined from the monopole via 
Eq. up to the ellipticity amplitude. We obtain an ap¬ 
proximate mass distribution for this cluster sample via the 
mass-richness relation of Rykoff et al. (2012J, 


M 200 = 10 14 M Q /h x exp (1.72 + 1.08 In (A/60.)), (24) 

where A is the cluster richness. We expect this relation to 
only be a rough approximation for our sample: it was derived 
for maxBCG clusters and involved assuming a cosmology in 
order to apply abundance matching techniques. The mean 
mass is M 200 ~ 3 x 10 14 M©//i, roughly an order of mag¬ 
nitude larger than the mean mass of our LRGs. To obtain 
the concentration we use the mass-concentration relation of 
Duffy et al. (2008), yielding c ~ 4.7 at the mean mass. Fi¬ 
nally, the mean redshift of the sample is z « 0.23. Note 
that many of these redMaPPer centrals are in fact LRGs 


(Hoshino et al. 2015): this measurement is therefore not com¬ 


pletely independent from our LRG measurements. LRGs oc¬ 
cupy a wide range of halo masses, and those which are also 
redMaPPer centrals will occupy some of the largest halos. 
Since the lensing quadrupole is expected to be proportional 
to the monopole, this measurement of cluster quadrupoles 
should have a higher signal than the LRG measurement: as 
will be shown below, the data matches this prediction. 

The result of the measurement and best-fit model is 
shown in Fig. [5] The constraints on the systematic shear 
permit a larger systematic amplitude than for LRGs, but 
the measurement is still signal-dominated. For example, at 
R rsj 0.1 Mpc /h the systematic amplitude is 20% (35%) of 


the halo ellipticity signal for the AE^ ’ (AE^ ') estimator. 
Reassuringly, the amplitudes of AE^ and AE^ are sim¬ 
ilar to each other, as expected if this shear component is 
insensitive to this systematic. 

The expected amplitude of the lensing signal is larger 
for this sample than the LRGs, due to the larger mean mass 
of clusters. In the absence of systematics such as misalign¬ 
ment of light and dark matter and correlation of lens and 
source ellipticities, we would expect the measured signal for 
clusters to be ~ 2 — 3 times larger, due to the 10 times larger 
mean mass. (Recall that Eq. <© requires the quadrupole to 
be proportional to the monopole, which we model with an 
NFW profile.) The observed difference is close to this, but 
the detailed fits of ellipticity show a slightly smaller value for 
the clusters. Although higher signal-to-noise measurements 
will be needed to tell whether this ellipticity difference is sig¬ 
nificant, it is plausible that these systematics are stronger for 
clusters. With our cut P cen > 90%, the redMaPPer galaxies 
should mostly be well-centered within the halo; nonetheless, 
it is still hard to compete with the LRG sample for which 
there is almost always only one very bright galaxy, making 
the halo center very easy to determine. Also, it may be that 


the cluster halo’s major axis correlates better with the ma¬ 
jor axis of the member galaxy distribution. However, this is 
a very different method susceptible to different systematics 
(Evans & Bridle 2009), and so we leave it for future work. 

Our fits of the systematic shear (top panel of Fig. |6|, 
show that this systematic is larger (1 — 2 a significance) for 
redMaPPer clusters than LRGs. The fits for halo ellipticity 
e are also shown in Fig. [6] The best-fit ellipticity is e = 
0.21, corresponding to an axis ratio 0.81. The reduced of 
the best-fit model is 20.6/18, an acceptable fit to the data. 
The hypothesis of no halo ellipticity is ruled out at slightly 
larger than the 3 a confidence level. The main reason for the 
smaller detection significance is that the number of clusters 
is less than 5% of the number of LRGs. As discussed above, 
the best-fit ellipticity may be smaller than LRGs due to 
greater misalignment of light and dark matter, but we leave 
a detailed study of major axis estimators for future work. 


5 SYSTEMATIC AND OTHER CHECKS 
5.1 Systematic Lens-Source Alignment 


Since the systematic described in § |2.2| is due to camera 
and PSF related effects, it should be present regardless of 
whether a given source galaxy is at higher redshift than the 
LRG lens. This suggests a model independent test of the 
systematic alignment between lenses and sources, in which 
we only use foreground galaxies as “sources.” Such galaxies 
will not be lensed and yet should be affected similarly by 
the systematic, providing a clean measurement of the sys¬ 
tematic alignment. We perform the measurement using all 
foreground sources with redshifts (z s ) < zl — 0.1, where (z s ) 
is the mean of the individual source’s redshift distribution 
P(z s ). Note that since the E cr i t weighting is irrelevant for 
sources with redshift below the lens redshift, we drop this 
weighting and measure 7 instead of AS. One qualification 
in interpreting this test is that changing the redshift cut 
on the “source” sample can also change other mean proper¬ 
ties such as size and surface brightness, which correlate with 
the level of PSF related systematics. Furthermore, high red¬ 
shift LRGs will be weighted more highly in this test than in 
the actual measurement, since they have more low redshift 
“sources” available. 

The results are shown in Fig. [7] Qualitatively it is con¬ 
sistent with our chosen systematic model in Eq. (16): the 


imaging systematics cause a positive AE^ ^ signal and 
do not contribute to AE^For a quantitative check we 
can scale this plot by the average E cr it factor of our mea¬ 
surements, E cr it ~ 8000 M 0 /i/pc 2 . The magnitude of the 
best-fit power law systematic in Fig. [7] is ~ 0.001 zb 0.0005 
at R rsj 0.1 Mpc /h. Then multiplying by 8000M©/i/pc 2 we 
have 8±4M 0 /i/pc 2 at R ~ 0.1 Mpc /h. Thus, the fit is only 
2cr away from zero, with the mean about four times the es¬ 
timate from the upper-edge of the la contour of Fig. [3] If 
we allow for negative amplitudes of the systematic, the level 
shifts lower to be consistent with zero at ~ 0.0007 d= 0.0007. 
In this case the mean systematic is between two and three 
times higher than the la edge of Fig. [3] While this is a higher 
value for the systematic, note that only the first 2-3 black 
points are significantly positive. The estimated systematic 
from this method could then be aggressively removed by 
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Figure 6. Same as Fig. [3] but for redMaPPer clusters. With a 
factor of 20 fewer objects, the noise of this cluster measurement 
is considerably larger, but the significance is still above 3cr due to 
the larger signal. 



Figure 7. Same as Fig. [2] but uses foreground sources which are 
not lensed to measure systematic alignment of LRGs and fore¬ 
ground galaxies. The best-fit power law systematic (solid line) is 
fit to the AEi points, while AE 2 is consistent with zero. 


discarding the first few points in the measurement. We have 
tested the effect of doing so, and find that discarding the 
first 2 (3) points for all four estimators results in a detection 
significance of 3.4cr (2.3cr). Since the detection significance 
with all points included is 4a, there does not appear to be 
any significant contamination from possible systematics at 
small R. 


5.2 Galaxy-Halo Misalignment 


Another potential systematic that could dilute the halo axis 
ratio measurements is misalignment between the position 
angles of the LRG light and dark matter distributions. While 
our upper bounds on axis ratios are robust even in the pres¬ 
ence of misalignment, the best-fit constraints are not. If our 
best-fit axis ratio of 0.78 is large compared to expectations 
from theory and other measurements, then misalignment is 
a possible cause. Adhikari et al. ( 2015| found axis ratios of 
q ~ 0.65 for halos with mass ~ 10 M©//i, falling to ~ 0.6 

for their largest mass bin of ~ 5 x 10 12 M©/h. Our best- 
fit axis ratio value is significantly larger than 0.6, but their 
most massive bin is still about a factor of 10 smaller than 
our LRG halos, so it is difficult to make a direct compari¬ 


son. In addition, Oguri et al. (2010) measured a best-fit axis 


ratio of 0.54 for 18 clusters specifically chosen to be well- 
fit by an elliptical NFW profile. Being very massive and 
specifically chosen to give good ellipticity constraints, these 
clusters may give a lower boundary on the expected halo 
axis ratios. Given the additional uncertainty in galaxy-halo 
alignment distributions, and the interplay of this systematic 
with axis ratio constraints, we leave a careful comparison 
to theory for future work. For now we only derive the ef¬ 
fect of galaxy-halo misalignment on our measurements and 
compare two possible alignment distributions. 

Any misalignment will lower the measured signal, thus 
requiring a smaller best-fit value of e for the model to com¬ 
pensate. If the angle between the galaxy and halo major 
axes is 0 O ff, then Eqs. (|4)[5| become 


Sent 7+ = (e/2) poryo ~ h~ h\ cos 2(0 - 0 off ) (25) 

E cri t 7x = (e/2) [-Ii + I 2 \ sin 2(0 - 0 off ). (26) 

Expanding the cos 2(0 — 0 o r) and sin 2(0 — 0 o ff) factors 
yields terms proportional to cos20 o ff and sin20 o ff. Insert- 
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10 Clampitt and Jain 


ing into Eqs. (|9[ |ll| ) and integrating over the distribution of 
alignments P a iign(0 o ff), the sin20 o ff terms average to zero. 
(Halos are equally likely to be misaligned on either side of 
the galaxy major axis, so that P a ii g n(0off) = Malign(—0 o ff) by 
symmetry.) This leaves 

(X cr it7i) = T>(e/4) [(2/i — Xo?7o) cos4^ + 2/2 — Xo?7o] (27) 
(X cri t 72 > « D{e/A) [2 h ~ X 0 ry 0 ] sin 4(9 , (28) 

where the dilution factor due to misalignment is given by 


5.3 Miscellaneous Tests 


D = 


nn /2 

/ d^offPalign(^off) COS 20 o ff . 

J — 7r/2 


(29) 


Thus misalignment is degenerate with halo ellipticity. 

In order to gauge the range of possibilities for misalign¬ 
ment, we compare various models to see how each would 
dilute the signal. The measurements of the large scale cor¬ 


relations of SDSS LRG’s by Okumura et al. (2009) show a 


signal over a wide range of scales (larger than those investi¬ 
gated here). By comparison with CDM dark matter halos in 
N-body simulations, the authors interpret the results as con¬ 
sistent with a misalignment of the LRG light with its parent 
halo, modeled as a Gaussian with a dispersion of a = 35°: 


Palign(0off) OC C off 


— 0 off / 2cr 


(30) 


The measurements of Li et al.| (2013) show that this broad 
Gaussian distribution also helps to explain the alignment of 
CMASS galaxies with large-scale-structure. Assuming this 
misalignment distribution, we find a large dilution factor 
D ~ 0.50. Applying this misalignment distribution of Oku¬ 
mura et al. (2009) to our measurement, we obtain a best-fit 


axis ratio of ~ 0.6, consistent with the simulations of Ad- 


hikari et al. (2015) 


However, a contrasting result is given by Ten neti et al.| 


(2014 1 ), which used hydrodynamical simulations to study the 


alignment of dark matter halos with their central galaxies. 
For their largest mass bin with Mhaio > 10 13 M©//i at z — 
0.3, the sample closest to our LRGs, Tenneti et al. (2014) 
found the 2D alignment distribution 


Palign(0off) OC e 9 - 28 l 6 ' off l 


(31) 


accurately describes their simulation results. This distribu¬ 
tion gives a much less drastic misalignment dilution, D « 
0.95. Note however, that the results between different hy¬ 
drodynamical codes have not converged: details of the sub¬ 
grid physics can still affect results substantially, including 


results for the alignment of galaxies and halos (Velliscig et al. 


2015). Nonetheless, if the above distributions span the range 


of possibilities, our best fit axis ratio would fall somewhere 
between 0.58 — 0.77 (ellipticity between e = 0.25 — 0.48). 

In the presence of this theoretical uncertainty, another 
approach is to attempt three-point lensing methods such as 


the estimator of Adhikari et al. (2015), which should produce 


constraints on ellipticity that are immune to galaxy-halo 


misalignment. (See also Simon et al. (2012).) However, their 
method requires much more data than is currently available: 
the authors find that LSST and Euclid should provide strong 
detections, but ongoing surveys such as PanSTARRS, DES, 
and HSC may have to wait until completion before having 
a chance at a detection of halo ellipticity with this method. 


We have attempted a few more measurement tests, which we 
summarize briefly here. First, as noted in § 2.1 Eqs. (|9| |11| 


suggest estimators that involve weighting AX2 by sin (4 6) 
and AXi by cos (40). We attempted to measure with such 
estimators but found no improvement in the results. Not 
surprisingly, constraints from AXi were degraded due to 
averaging out the terms without 6 dependence in Eq. 
There are no such terms in Eq. (11); however the elliptic¬ 


ity constraints from AX2 shifted downward by la compared 
to our fiducial constraints. The choice of weights and esti¬ 
mators merits further study. Second, our AX^ and AX^ 
estimators will have some contamination from cosmic shear, 
which aligns nearby sources (since £+ >0). This means cos¬ 
mic shear acts with opposite sign to our elliptical halo lens¬ 
ing and with the same sign as the systematic shear that 
we have treated in detail in Secs. lT2l and 01 Since cosmic 
shear also has roughly the 1 /R dependence of our systematic 
model (see for example Fig. 8 of |Lin et al.| (2012)), to some 
extent it is already folded into our systematics treatment. 
Third, the lens galaxies for which we have shape measure¬ 
ments generally have very small shear measurement errors 
and a significantly non-zero ellipticity magnitude. We tried 
cutting on both these quantities, removing the few lenses 
that are poorly measured or very round. However, neither 
cut affected the results significantly, likely since only a small 
fraction of lenses were removed by these cuts. This issue may 
be more important for the full galaxy population, which in¬ 
cludes much fainter lens galaxies. Finally, we tested the im¬ 
pact of satellite LRGs since they may have poorer alignment 
of light with halo mass. We found that removing the ~10 
percent of LRGs in multiple systems did not change our re¬ 
sults (the best fit ellipticity value shifted upwards by about 
half-sigma). 


6 CONCLUSIONS 

Using a large sample of SDSS LRGs, we have studied the 
quadrupole shear of the dark matter halos and ruled out 
the hypothesis of no ellipticity at the 4 a confidence level, 
and ellipticity < 0.12 at the 2 a level. Our observable uses 
carefully chosen angular averages of the two components of 
Cartesian shear, none of which have a contribution from 
the usual circularly averaged monopole term. Our estimator 
and modeling have advantages over previous results with 
tangential- and cross-shear estimators. One component of 
our Cartesian estimator is insensitive to the primary sys¬ 
tematic in the literature ( Mandelbaum et al.|2 006), a spuri¬ 
ous alignment of lens and source ellipticities, allowing us to 
isolate the shear from the halo quadrupole. We obtain un¬ 
biased measurements from the component that is affected 
by the systematic by modeling and marginalizing over its 
effect, eliminating the need for subtracting it off and hence 
increasing noise. 

Our key result on the halo ellipticity of LRG’s is pre¬ 
sented in Fig. [3] The central value of the LRG halo ellipitc- 
ity is 0.24 =L 0.06 corresponding to an axis ratio of 0.78. We 
find that our measurements agree with simulation studies 
of CDM halos provided one folds in a level of misalignment 
of mass and light. A Gaussian distribution of misalignment 
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angles with a — 35 degrees, suggested by |Okumura et al.| 
(2009), fits our results (though is higher than the recent 
simulation-based results of |Tenneti et al. (2014)). We also 
measure the ellipticy of redMaPPer clusters and find results 
consistent with those of LRG halos, though with slightly 
smaller statistical significance owing to the smaller sample 
of suitable clusters. 

Our modeling could be further improved by simul¬ 
taneously fitting for halo mass and ellipticity and then 
marginalizing over halo mass. The current model has some 
uncertainty since we use a fixed mass and the predicted 
quadrupole signal is proportional to the surface density So, 
which is larger for more massive halos. The la mass uncer¬ 
tainty in Mandelbaum et al. (2008a) for the LRG sample is 
~ 10%. Consider the case where the true halo mass is 10% 
higher than our assumed mass of 4 x 10 13 M©//i. In that 
case we find a true best-fit ellipticity e = 0.22: thus under 
these assumptions our fiducial best-fit value of e = 0.24 is 
an overestimate of about 8%. 

While marginalizing over mass would make us less sus¬ 
ceptible to such biases, it would likely degrade our con¬ 
straints somewhat relative to our current model with a sin¬ 
gle fixed mass and concentration. This is because the same 
background sources would be used to determine both the 
monopole and quadrupole parameters, and the shape noise 
would then be slightly correlated between quadrupole and 
monopole elements of our combined data vector. However, 
since the monopole signal is nulled with our quadrupole es¬ 
timator, the signal itself would add no additional covari¬ 
ance. The dependence of the signal on concentration is much 
weaker, and so marginalization over it is likely to be less im¬ 
portant. Nonetheless, concentration may be worth marginal¬ 
izing over along with higher order parameters such as mass- 
concentration scatter and mass-luminosity scatter in future, 
higher signal-to-noise measurements. 

Ongoing surveys (DES, KiDS, HSC) have superior im¬ 
age quality and depth compared to SDSS. With the poten¬ 
tial for significantly smaller errors on halo ellipticity mea¬ 
surements, it will also be necessary to explore additional 
tests for systematics and improve the modeling compared to 
the analysis from SDSS presented here. In the previous sec¬ 
tion we have summarized a set of tests that can be pursued 
in greater detail. Other extensions could include separating 
central lens galaxies; using other optical estimators of the 
alignment of the halo such as the distribution of neighboring 
galaxies; exploring other estimators for the halo ellipticity; 
and extending the measurements to field galaxies to study 
the trends with luminosity and galaxy type. Also in the fu¬ 


ture, higher-order lensing effects such as flexion ( Hawken fc 
Bridle 2009) will help to increase the signal-to-noise of halo 


ellipticity measurements. 


6.1 Comparison to other work 

In §|4.3|we com pared our results in some depth with those of 


Mandelbaum et al. (2006). Now we briefly survey the other 


attemp ts to measure halo ellipticity from data. |Parker et al.| 
(2007) used a sample of ~ 2 x 10 5 lenses and ~ 1.3 x 10 6 
sources in 22 sq. deg. of single-band CFHTLS data to at¬ 
tempt a halo ellipticity measurement. Without photomet¬ 
ric redshifts the split between lenses and sources was done 
based simply on galaxy magnitude in the one available band. 



x [Mpc//i] 


Figure 8. The ellipticity of the dark matter distribution (black 
ellipses) is shown over the range of scales sensitive to our mea¬ 
surement of quadrupole lensing. The visible baryons (red ellipse 
at center) are essentially confined to scales ~ 40 times smaller. 
Regardless of the details of the theory, any attempt to modify 
gravity and use just the baryons to describe lensing will produce 
rounder contours with increasing radius. 


They found a 2a hint of elliptical halos by looking at the ra¬ 
tio of tangential shear along the major and minor axes of 
the lenses. With data from the second Red-sequence Clus¬ 
ter Survey (RCS2) |van Uitert et al. (2012) used an estimator 
involving the ratio of tangential shear along the lens’ major 
and minor axes, as well as the tangential- and cross-shear 
estimators Eqs. (22| 23), and none gave a significant detec¬ 
tion. (See also Hoekstra et al. (2004) for an earlier study 
of RCS.) There has been more success in the cluster lens¬ 
ing regime, where the intrinsic signal of elliptical halos is 
higher. Using 4,300 cluster lenses in SDSS |Evans &"T5 ridle 
(2009) obtained a 2.7cr detection of signal. Furthermore, us¬ 
ing a subset of 18 clusters that are well fit individually by 
an elliptical halo model, Oguri et al. (2010) estimated a 7 a 
detection of ellipticity for those specific clusters. 


6.2 Dark matter and gravity theories 

As we have seen, galaxy-galaxy lensing provides information 
on halo shapes over a wide range of distances. It turns out 
that this range is very useful for constraining some theories 
of dark matter and gravity. Typically the virial radii of these 
galaxies are at least ten times larger than their visible sizes. 
Lensing measurements of the kind we have presented span 
this range of distances, in the case of our LRG sample it 
covers the range from 0.1 Mpc /h to 1 Mpc /h and beyond 
with the best signal-to-noise at about 0.5 Mpc /h. We discuss 
next two examples of theories that can be tested by our 
measurements and future improvements. 

The lower bounds on halo ellipticities can constrain the¬ 
ories that attempt to do away with particle dark matter by 
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modifying gravity. In Fig. [8] we show the ellipticity of the 
dark matter distribution from 0.1 Mpc /h to 2 Mpc /h, the 
range of scales sensitive to our measurement of quadrupole 
lensing. The visible baryons, also pictured, are essentially 
confined within 0.05 Mpc /h. Regardless of the details of 
the theory, any attempt to modify gravity and use just the 
baryons to describe lensing will produce rounder contours 
with increasing radius. Hence for MOND-like gravity the¬ 
ories, including recent extensions like |Khoury (2015), our 
measurements pose a challenge. In such theories, the ex¬ 
ternal tidal field can produce a quadrupole dependence on 
large scales, but it would be very difficult to obtain the right 
scale dependence over the range of our measurements. On 
the other hand, halos in GR-based CDM cosmologies have a 
roughly scale-independent ellipticity that is consistent with 
our measurements. 

Another theoretical scenario of interest is self¬ 
interacting dark matter. If the dark matter particle has in¬ 
teractions that lead to scatterings with order unity proba¬ 
bility over the age of the halo, then these scatterings make 
the halo rounder. Since the density is higher in the central 
regions, the effect is stronger closer to the center of the halo 
( Peter et al.|2013 ). Thus there is a scale dependence that dis¬ 
tinguishes self-interacting dark matter halos from standard 
CDM ones. To the extent that our ellipticity measurements 
match CDM predictions, we can constrain the interaction 
cross-section for self-interacting scenarios. Figure 3 of [Peter] 
|et al. (2013|) shows the variation of the axis ratio over 0.1 to 
1 Mpc/h, a range well probed by our measurements. See also 
the recent study of Kaplinghat et al. (2014) which examines 
how baryons, which may dominate in the inner parts of the 
galaxy, can counter the isotropizing effect of dark matter 
self-interactions. 

We leave the connection of our measurements to these 
theoretical scenarios for future work since it requires a care¬ 
fully study of the scale dependence of halo ellipticity. It may 
be that we will need larger datasets with better lensing sig¬ 
nal as expected from ongoing surveys to obtain meaningful 
constraints. 
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